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The study of the effect of resonant expansion systems on the 
propagation of long waves is presented in this thesis as a design 
problem in coastal hydraulics. Two previous studies which have been 
made on similar systems are presented in Chapter I to provide back- 
ground material on the subject. The first study presented was made 
by Lamb (1916); and it demonstrates that reflection of long waves can 
be obtained by an expansion system and that the amount of wave energy 
that is reflected is a function of the geometry of the system. The 
second study presented was made by Valembois (1953) and is based 
on a hydrodynamic impedance theory. An important point to note in 
Chapter I is that both of these studies use a scalar equation of pressure 
continuity and a vector equation for the conservation of mass to obtain 
a solution. A new solution is derived in Chapter 2 by formulating a 
boundary value problem which incorporates the same type of boundary 
condition equations noted earlier. The problem is formulated in con- 
siderable detail to demonstrate the type of problem solution technique 
which is required to solve problems arising in coastal hydraulics. The 
details given in the problem statement become important later when 
the effect of geometrical changes of the resonant system are evaluated. 


Chapter 3 discusses the interdependance required between Civil 











Engineering and Oceanography to effect a complete solution to a design 
problem in coastal hydraulics. Several methods for presenting wave 
data are given, and one of these methods is selected for the design 
problem being considered. Finally, a linear analysis is employed in 
Chapter 4 to evaluate the effect of the resonator on the propagation of 
long waves by combining the solution derived in Chapter 2 with the 
spectral method of wave data presentation selected in Chapter 3. 

The results of the linear analysis are presented in graphical form 

as a measure of the amount of wave energy that is reflected and 


transmitted by the resonator. 
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The response of a single resonant expansion chambertomthe 
periodic pressure fluctuations of ocean waves is to act as a rigid 
vertical reflecting surface. The effective frequency band-width over 
which reflection is significant may be determined for a given harbor 
geometry by means of a graphical linear analysis incorporating the 
reflection or transmission coefficient and wave data presented by a 
Bretschneider Power Sprectrum. The most critical dimension of 
the resonator was determined to be the length of the expansion, La: 
measured transverse to the centerline of the channel; the effect of 
the resonator was maximized when this dimension was either one- 
quarter or three-quarters of the design wavelength, L,. The influ- 
ence of the width of the resonator, Dp: measured parallel to the 
centerline of the channel was determined to be maximal when this 
dimension was one-half of the design wavelength, L,,. ihemetiec timc - 
ness of the resonator is also a function of the width of the main 
channel, VE and this parameter was incorporated in the expressions 
for the coefficients. 

Because of its analogy to electrical band-stop filters, the 
possibility for improving the effect of the resonator by constructing 
several resonators in cascade appears good and offers an opportunity 


for extending the results obtained here. 
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As the logistical requirements of nations grow and expand, new 
techniques for designing harbors will be required to accomodate the 
expansion of littoral ports of communications. The port areas which 
are sheltered naturally from the ocean forces and which require little 
or no additional protective structures have already been exploited. 
Designers, therefore, must address themselves to the problems of 
the less protected areas where this expansion must occur. One of 
the major concerns in an unsheltered area is the design of breakwaters. 

The purpose of breakwaters is to cause a reduction of wave 
heights in its lee (Wiegel, 1964). The conventional methods used for 
sheltering the saieanee of a harbor between breakwaters may be 
divided into three peneralmmeecs(Miourec. |b) a altoree methods 
shown require shipping vessels to maneuver in order to enter or to 
leave the harbor. This required maneuvering is both expensive and 
inconvenient to the user. In addition, the third method shown requires 
a deep shipping channel close to the shore, and may require the addi- 
tional expense of periodic dredging. The alternatives which have been 
considered to provide a direct access to the shipping channels from 
the harbor have included pneumatic and hydraulic devices as well as 
resonant expansion chambers (Wiegel, 1964). The development of the 


latter is considered here. 
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Figure |.l. Three conventional methods used to protect 
breakwater entrances. 
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Wave Reflection by an Abrupt Expansion 


Lamb (1916) investigated the partial reflection of a wave 
resulting from an abrupt change in the cross section of a channel. 
His results were based ona linear solution to the wave equation. 

(A linear solution is presented in Chapter 2). 

A right-handed cartesian coordinate system (z-axis positive 
upwards) is oriented along the centerline of a channel at the mean- 
water surface with the origin located at the discontinuity (Figure 1. 2). 
The wave amplitudes and relative velocities to the left of the discon- 


tinuity are given by 


=e ee (1.1) 
1 1 
ae 2° 3 (ee ye ea (1. 2) 
mel ai | 


and to the right of the discontinuity by 


Ny = 9(t- aaa es 
2 

— _ 2 .. . 2S: 

iy: = a p (t c,! ’ ee) 


where F is an arbitrary periodic function representing the incident 
wave; f and @ are arbitrary periodic functions representing the 
reflected and transmitted wave, respectively; and cy and c, are 


the wave celerities of the wave systems to the left and to the right of 


the discontinuity, respectively. The linear superposition theory is 
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Figure 1.2. Partial wave reflection by an abrupt expansion 
in ag channel. 








valid provided that the small amplitude assumption holds; i.e., 





<< |  hOa = leer (1. 5) 


The assumption that the waves are propagating in shallow water 


requires that 





See (lee6') 


where L, is the wavelength in deepwater (Kinsman, 1965). By invoking 
the shallow water wave assumption given in Eq. (1.6), the wave 


celerities may be approximated by 
c = gh. tor 1 = lez. (IF, °7) 


There are two boundary conditions which must be satisfied at the 


Origin (x=y=0). The conservation of mass requires that 


Bh ih 


"| 1 See 2 oe oe 


where Bi: B, are the breadths of the channel measured at the mean- 


water surface, and h_., h. are the mean-water depths. The assumption 


l Z 
that the fluid motion will be sensibly uniform along and parallel to the 
centerline of the channel for small distances (comparedtoa wavelength) on 
either side of the discontinuity neatnre 2 that there be no sensible 
change in the mean-water surface across the discontinuity. This con- 


dition requires that the continuity of pressure across the abrupt change 


be 








Le al (1.99 


Substituting Eqs. (1: 2) and (124) into es) scives 


Bh 
(1. 10) 


Bh 
2 
eo (1s) Mrs Weel Oe 
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] 
and the substitution of Eqs. (1.1) and (1.3) into (1.9) yields 


(ert } 


F (t) + £(t)= o(t) at x= 0. 


The ratios of the reflected wave and of the transmitted wave to the 


incident wave determine the reflection and transmission coefficients, 


respectively; iwe., 


j _p 
KrF and K= 


These coefficients are determined by the simultaneous solution of Eqs. 


(1.10) and (1.11). By means of Eq. (1.7), the reflection and trans- 


mission coefficients, respectively, are 


Paes nem — 
a af ¢ 
R Bic, Boc., 
and 
K_= mete (leak 3) 
ae Bic) Bc, 


The energy contained in the reflected and in the transmitted waves is 


equal to the energy of the incident wave if the sum of the squares of the 


two coefficients is equal to unity; i.e., 
(1. 14) 


2 2 
K 2s ip 
Rt Be 


Equations (1.12) and (1.13) establish that the effectiveness of the 








abrupt expansion for reflecting waves is dependent on the geometry of 
the channel. 

A similar procedure to the above approximation given by Lamb 
(1916) could be employed to evaluate the reflection and transmission 
resulting from an abrupt contraction of a channel and the two solutions 
superimposed in a limiting case. However, the basic assumption that 
uniform conditions exist parallel to the centerline of the channel would 
now introduce a greater error of approximation, as the dimension of 
the length of the expansion measured parallel to the centerline of the 
channel approached small multiples of the incident wavelength. These 
results obtained by Lamb (1916) indicate that substantial reflection 
occurs from an abrupt expansion in a channel and that a further investi- 


gation of an expansion system is justified. 
The Resonant Expansion Chamber 


The resonant expansion chamber (resonator) takes advantage 
of the periodic character of the surface disturbances and of the pres- 
sure fluctuations produced by a wave (Valembois, 1953). The solution 
derived by Valembois includes a wave system transverse to the center- 
line of the channel which is induced in the resonator; this wave system 
was neglected by Lamb in the preceding example. The derivation again 
assumes a linear solution to the wave equation but provides an additional 


parameter for evaluating any energy dissipation which may result. The 











measure of energy dissipation utilizes the Neumann impedance theory 

for a hydrodynamic oscillatory system; this theory is analogous to 

the impedance theory for electrical oscillatory systems (Defant, 1961). 
The hydraulic impedance of the discontinuity (Figure 1.3) is 

in accordance with Ohm's law for the analogous electrical oscillatory 

system. The effect of the resonator is to alternately withdraw from 

and discharge into the main channel a quantity of fluid which is pro- 

portional to the surface distortion associated with the incident wave 

system. The constant of proportionality is a measure of the hydraulic 

impedance of the resonator; i.e., 


A 


yh = 
V 





ls Ve 


A right-handed cartesian coordinate system (z-axis positive upwards) 
is Oriented at the orthogonal intersection of the centerlines of the main 
channel and of the resonator at the mean-water depth. The linear 
shallow water approximation for the wave celerity may be written as 
in the preceding example as 

Bo Nila. (1. 16) 
where h is the mean-water depth of the channel. The relative velocity 
U(x, t) and superelevation H(x,t) in the section to the left of the dis- 


continuity at A are given by the following relationships; 


U Ge) - Re eae) (1.17) 
Cc Cc 


NG) erat © (etn ede (1. 18) 
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Figure 1.3. Action of a resonator on the propagation of a 
gravity wave. (after Valembois, I953) 
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and to the right of A by the following: 


U’ (x,t) = F° = ) aie) 
Eis ie Bolte = (1. 20) 


where F is an arbitrary periodic function representing the incident 
wave; and f,F* are arbitrary periodic functions representing the 
reflected and transmitted waves, respectively. Since the medium of 
propagation is assumed to be at rest, all relative velocities will denote 
absolute velocities. The channel is also assumed to be infinite in 
extent, both in the upstream and downstream directions; therefore, 

no reflection of the reflected wave or of the transmitted wave by a 
terminal boundary is considered. 

The desired reflection and transmission coefficients will be 
determined exactly as before by considering the simultaneous solution 
of two boundary conditions at A: the continuity of pressure and the 
continuity of flow. The first boundary condition requires that 

ss ee Rt (1. 21) 


A A’ 


and the latter requires that 


n A a a (lee 2 } 
pubstituting as. (1218) and { lLo20) into 1.21) yields 


me pf Re aoe 


eez 
A A A Cc A’ ( >) 


and substituting Eqs. (1.17) and (1.19) into (1. 22) yields 














Al 


ee ea | (1) 2) 


Substituting Eq. (1915) foreveand designating the teculiing ratio by 


of EES Ge | (1. 25) 


the simultaneous solutions of Eqs. (1. 23) and (1. 24) result in the 


following values for the reflection and transmission coefficients, 


respectively: 
K_ = | . | (1. 26) 
Ig fe 08 oa ah 
and 
(0) 
Ke Share ae falipneZ 7) 


The sum of the squares of the two coefficients is not equal to unity 

as in the previous example by: Lamb (1916), except when the value of 
@ (and hence Z) is exactly equal to zero. The requirement that 
shewalwexof Ze giawensby Eq. (&. | 5)ebesequal toszeroscommespondssto 
no surface elevation in the resonator (e. g the resonator discharges 
into a large reservoir) or to an infinitely high velocity of discharge. 
The hydraulic impedance may be written in a form analogous to the 
electrical impedance; me. , 


1 


A= rer jee UO au 





i (1. 28) 


where r is proportional to the energy dissipative forces; and £, y 
are proportional to the inertial and potential (e.g., gravity) forces. 


Requiring Z to vanish is analogous to a tuned resonant circuit in which 











12 


i (1. 29) 


where O is the angular frequency of oscillation. The value of Z will 
vanish at resonance only in an inviscid fluid (or in an analogous con- 
servative system). Equations (1. 26) and (1.27) may be used for 
evaluating the amount of energy dissipation and the corresponding 
amount of reflection of waves propagating ina real fluid only whena 
value for r is known or may be determined. 

The difficulties enecunuaeee in the application of the results 
of Valembois (1953) lie in the evaluation of the parameter r and in 
the criteria for dimensioning a tuned resonator. The values of the 
hydraulic resistance and reactances are very small and are usually 
determined empirically. Secondly, an important omission in the 
theory is a rigorous analytical procedure for dimensioning the geo- 
Igetry of the resonator. The only reference given to a dimension was 
the empirically determined value for the length of the lateral expansion 
which was equal to one-quarter of a wavelength. Since surface gravity 
waves have no natural frequency of oscillation to which the natural 
frequency of the resonator may be matched, the electrical analogy 
cannot be employed as a means for dimensioning the resonator. 

This thesis is addressed to the evaluation of a reflection and 
of a transmission coefficient which are independent of any empirical 
data required for measuring the small dissipative forces and which 


relate any geometric changes in the resonator to its effectiveness for 





Ve} 


reflecting waves in terms Of thesebaracteristic parameter olor 


incident wavelength. 











éve THE SOL LION 


Formulation of the Boundary Value Problem 


The derivation of the desired coefficients will proceed ina 
manner similar to that of the examples presented in Chapter 1 by the 
formulation of a boundary-value problem. A right-handed cartesian 
coordinate system (z-axis positive upwards) is oriented at the orthogo- 
nal intersection of the centerlines of the channel and of the resonator 
at the mean-water level (Figure 2.1). The linear second-order homo- 
geneous hyperbolic partial differential equation which approximately 
describes the instantaneous surface elevation of the fluid in the domain 
mutually occupied by the channel and by the resonator is the wave 


equation; i.e., 





—" =< an. (2° ap 


where k has a range of values from unity to three. The wave form 
propagates at a constant celerity of 


© a (Z.,a) 


Neglecting the non-linear effects of the boundary conditions requires 
that the following assumptions be made (McLellan, 1965): 
(1) The amplitude of the surface disturbance is 


small compared both to the wavelength and to 
the depth of the fluid. 
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Figure 2.1. The wave systems considered in evaluating the 
effect of the resonator on long waves. 
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(2) The channel is of uniform depth. 


(3) The fluid is inviscid and irrotational. 

( 4) The fluid is incompressible and homogeneous. 

(5) Coriolis acceleration may be neglected. 

(6) Surface tension may be neglected. 

(7) The rigid boundaries are smooth and impermeable. 
(8) The atmospheric pressure along the water surface 


is constant and uniform. 
The solution to Eq.( 2.1) for one dimensional space may be 


expressed in functional notation by 


‘tele taas (Zug ds) 
The most general solution is given by D'Alembert's Theory and is of 


the form 


n (x t)8= e(les=eOq) ee Eutlae+ Oe (2. 4) 
mnere F is an arbitrary periodic surface disturbance propagating in 
the positive x direction and G is anarbitrary periodic surface dis- 
turbance propagating in the negative x direction. One simple-harmonic 


function which satisfies Eq. (2.4) is a cosine function; i.e., 


(2, +) = AL cos(k aS + Ot te) . (2.5) 


where k has a range from unity to five corresponding to the com- 
ponents shown in Figure 2.1 and j has a range from unity to two. 


The numerical indices for k are now replaced by alphabetical 
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indices shown in Figure 2.1 in the following manner: 


1 = Mi 
2. = Wi 
3 = Mt 
4= Ri 

5 Sai 


where the upper case letters M, R correspond to the main channel 
and to the resonator components, respectively; and the lower case 
letters i, r, t Gornespond)to thevinei dentyanemlected mand transmitted 
wave components, respectively. The ze components are replaced 

by x and y for the coordinate system shown in Figure 2.1. Equation 
(2.5) is a set of five wave component equations containing ten unknown 


A 


constants: viz., the five values for the wave amplitudes A, _., : 
Mi Mr 


A. _, A_., A and their associated phase angles 


Mt’ Ri’? “Rr Pe Mr eee 


a The set of equations expressed in Eq. (2. 5) are the following: 


Wi = Ay cos (Kx - Ot teQ,) (2. 6:1) 
Tur Aye COS (RX + Ott exe) (2. 6: 2) 
ce Yams OSS 5 alias) (2. 6: 3) 
Tei = An; cos(ky- oe + Ep;) (2. 6% 4) 
1 op Ont yt Ot + oe (2. 63 5) 


In a well-stated boundard-value problem, the unknowns contained 


in Eqs. (2. 6:1) may be evaluated by a given set of boundary and initial 








18 
conditions (Tikhonov and Samarskii, 1963). In the stating of the con- 
ditions required by the above to comply with the criteria for a well- 
stated boundary-value problem, the following assumptions are made: 


(9) The fluid medium of propagation is at rest 
(i.e., all velocities are absolute). 


(10) The terminal boundaries of the centerline of 
the main channel are sufficiently removed 
from the proximity of the resonator and, there- 
fore, are not required boundary conditions. 


{ i) Steady state conditions exist. 

(1:2) All geometric dimensions considered are 
small compared to wavelengths (i.e., long 
waves). 

C153} The cross sectional area of the channel and 


of the resonator consists of a horizontal 
bottom and parallel, vertical sides. 


(14) Shallow water conditions exist. 
Assumption (11) negates the requirement that the initial conditions 
be stated. Assumption (14) permits the use of the shallow water 
approximation for the wave celerity; i.e., 


Z 
c = gh. (Ze) 


The relative velocities shown in Figure 2.1 are related to Eqs. (2. 6:1) 


in the following manner; 








oC 
ee Mi M 
Mi . cos {kx - Ot+ VEL {Zee ) 
M 
A Cc 
= Mir Ma . : 
Oty = a cos (kx + Otte...) (2182.2) 











WY, 


A 
ae Mt M 
Ua a cos (kx - | (Jaca) 
M 
A c 
— R teeR . 
U Ri ee. ay cos (ky Otte.) (c.10s 4) 
R 
A res 
ea of) — . 
Ur, = = cos(ky + Otte, ) q (25.5) 


The boundary conditions will be prescribed along the free surface in 
terms of Eqs. (2.6:i) and along the vertical boundaries in terms of 
Mas. (2. 8:i). 

The free surface boundary condition prescribes the instantaneous 


continuity of pressure at the origin for all time: 


von) \inaelieeee gees Wars 


The vertical boundary conditions require two statements. The first is 


that the continuity of flow into tne domain is expressed by 


h i Ul h_= 0; 
ae [Ug + UE bphp= 9; (2. 0) 


—p> 


‘ <b , i 
own ayy in aa 


and the second is that the no flow condition across the rigid terminal 
boundary of the resonator is given by 
Uae 0. = 0. > ig 
pit Ved Oar i 
The above three boundary conditions are evaluated at the following 


coordinates: 


Eqert2. 9) at the onl aie (je. = = y = 0}: 


b 
Eq. (2.10) atx=a—, ya se; 


and Eq. (2.11) aty =-£., x= 0. 


R’ 
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Using the Boundary Conditions to Obtain a Solution 


The main channel incident wave may be arbitrarily selected as 
the reference wave profile by assumption (11) without any loss in 
generality; therefore, its associated phase angle is defined as 


=s0. sl 
Mi 0 (Zale) 


Eight independent boundary equations will result by the substitution of 
Eqs. (2. 6:1) and (2. 8i) into the three boundary conditions prescribed 
in Eqs. (2.9), (2.10) and (2.11) when the expressions obtained after 
substitution are expanded by means of the trignometric identities for 
the sine and cosine of the sums and differences of angles. The time 
dependence of the boundary conditions may be eliminated by equating 
the sine and cosine terms which result from the above expansion. The 


ten equations which result from the preceding algebra are the following: 


Ati + Air cosé,,. = Anat COS €,), (zee) 32.1) 
= Aatr sin Sow ee Anat sin ent (Z: 132.2) 
Ana aE Ve eos. ate Ar: cos €,; ck vee Bosh. .. (2 133) 
~ AN sin eur = tAR; sin epi Ans. sin eRy (2. 135-4) 
Antt 22S An; COS Ep, ~ An, COS Ep (25 luge 5) 
Anat sin eel Ap; sin Ri” Ap. sin eRe (2 saledinb ) 
Ani cos (kin- Ep) ~ fy cos ( kh, = ép,) (2e45..7) 








ZA 


- = = eS 8 | 
Ani sin (ke, Epi) Ay. sin (Ke, Ep) (25s 8) 
b, bp b 
——)-A k—— - -A k— aulieg 
Cy Oy Aang COSMK 5 ) mre?! 5 ne mt cos! a va (2 lee a) 
by b 
= - a +A k—- 
Sg es Tee gS CON gr) 
b b b 
R R R 
-G by HAL, Sink >") +A, Sink> Evie tA, ,Sin( k= 5 yin CZ lsu) 
BM b 
= CpPR LAR; sin{k =— - Ep;) a An. sin(k >— - Ep, ] 
byeimeans of Eqs. (Z.13:7) and(2.13: 8), it may be shown that 
= .14 
i RE iGo 
and‘that 
leis = + ; Le 
+2 n> ‘ri’ ‘re (2.2195) 


By squaring and adding Eqs. (2.13: 5) and(2. 13:6), Eq. (2.14) may 
be rewritten in the form 


A 
Mt 


Tey [2 cos( k£_)| * ema) 
R 
Equation (2.7) may be employed to equate the two celerities; i.e., 


a ia aha 7) 


Following the substitution of the known relationships and applying some 
additional trignometric identities, the reflection and transmission 


coefficients may be obtained by squaring and adding Eqs. (2. 13:9) and 








Zig 


(2.13:10). The reflection coefficient may be expressed by 





ey aN an 
sin (k>") + - )[sin(k-> ) - [tan (KZ, )| cos (k~> )] 
Phe Gen b b b b b b z 
R 2 
R R M M Re, M M2 
8 nce en S- pong eos + Ge) [sin(k7~ )-| tan(kg,)| cos (Kk “) )} | 
(2.18) 
and the transmission coefficient by 
UB 
7 cos ( , ) 
b b b By ale 


T b 2 ‘ 
M R 2 M M 
tS sin fin tan (Ke. ) |cos(k~")-]+ (>) [sin(k-> )-|tan(kg., )|cos(k>~)] 


as) 


The sum of the squares of Eqs. (2.18) and (2.19) is equal to 


unity and verifies that the energy in the system is conserved. 


x 


Evaluating the Geometric Parameters 


The criteria for dimensioning the resonator are the maximizing 
and minimizing of Eqs. (2.18) and (2.19), respectively. This minimizing 
and maximizing process requires the evaluation of the arguments of the 


teiemMoOmetric functions involving ¢_; b and Lev in terms of the tuned . 


R ke 


wave length, L,, and wave period, T,, of the resonator. 

The dimension which has the most critical influence on the 
value of the coefficients is the transverse length of the resonator, Lae 
which appears in the argument for the tangent function. The tangent 


function becomes discontinuous when its argument approaches a value 
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Oigaleaen S| a= radians, where n is any positive integer. The argu- 
ment reaches this value when 
A er 


i al 2 Uae la) ee tc. 20) 


where n must be equal to unity or two by assumption (12). Eq. (2. 20) 
agrees with the results obtained by Valembois (1953) when the value 
of one-half of the width of the main channel, vt is either a very small 
Bem cent of L, or equal to one-half L,,. The argument containing 

Q may be expressed in terms of the angular frequency, 0, by means 


R 


of the wave celerity equation; i.e., 








o 
ki, Se fp: (Zeaal) 
1s, 
Decreasing the dimension Ip in Eq. (2.21) toa value less than 7 
increases the tuned resonant frequency, ar » by a corresponding 
Lips 
factor; while increasing 1, to a value greater than = decreases 





the tuned resonant frequency, O , also by a corresponding factor; i.e., 
O 





a. 
_— > 
le i. < a then O 0 
JL 
aa 
and i R 7 , then O< 06° 


a 


where O is the tuned resonant frequency for Metimuw reflection. It 
may be seen from Eq. (2.16) that the resonator wave amplitude would 
become infinite at resonance for any finite value of the transmitted 
wave height, A,,... Since an infinite wave amplitude does not 


Mt 


- seem to be physically valid, the transmitted wave height, Aart 
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must be identicallyezero. 

The influence of the width of the resonator, Dp on the value 
of the coefficients is most easily evaluated by means of Eq. (2.19). 
Minimizing the transmission coefficient requires that the value of the 
cosine function in the numerator approach zero. The cosine function 
approaches zero when its argument approaches (2n ae , where n 
ls any positive integer. Therefore, the value of the width of the 
resonator must be 


ale 


pee -1) —— ; 
pote(2neeh) Fa 222) 


where n must be unity by assumption (12). If assumption (12) is 

valid, the value of b, determined bye hq. (2-22) 1S alrniasemamund: 

ence, bp may only decrease from this maximum value. The argu- 
ment containing DR may be expressed in terms of the angular frequency 


as before by 


b b 
R O R 
k ——_-_- =- —_ -—. eis) 
2 Cc Z (2523) 


Decreasing DR from its maximum value giveneby Eq. (2.22) tmereases 
the tuned resonant frequency We 

Finally,’a critical value for the width of the main channel, be 
is not easily evaluated by means of the equations for the coefficients; 
however, because the preceding derivation was based on the assumption 


that all dimensions were small compared to a wavelength, it may be 


determined that the width of the main channel must be a small per cent 











Z5 
of the tuned wavelength, Li. 

The evaluation of the geometry of the resonator system has 
been discussed by considering the effect of a geometrical change on 
the ability of the resonator to reflect an incident wave of a single fre- 
quency. Assumption (11) does not allow the application of the coeffi- 
cients to a constantly varying input inducing a transitory response. 
However, the experimental results obtained by Valembois (1953) 
indicate that appreciable reflection does occur over a range of fre- 
quencies. In fact, the response curves obtained by Valembois are 
similar to the response curves for electrical filter circuits (Olson, 
1958). The analogy of these two oscillatory systems introduces the 
possibility of determining an effective frequency band-width for which 


metlection may be significant. 








3.° PRESENTATION OF WAVE DATA 


The derivation of the reflection and transmission coefficients 
in Chapter 2 and their subsequent discussions dealt with an incident 
wave of a single frequency and wavelength. An analysis of an incident 
wave record reveals that waves in nature are composed of an infinite 
and continuous number of frequencies resulting from a complex 
phenomenon. Since the resonant structure is static, some means 
of determining a single frequency which is a characteristic measure 
of the complicated group phenomenon contained in the waves must be 
made as well as some means for evaluating the effect of the resonator 
when subjected to an input composed of an infinite series of incident 
frequencies. Fortunately, the solution to the first problem provides 
a means for resolving the second. 

The problems encountered in recording and in analyzing surface 
waves are currently under research by the discipline of Oceanography. 
A complete or even adequate treatment of these complex problems is 
beyond the scope of this thesis. The design and construction of the 
resonant breakwater serves to illustrate the interdependance required 
between Civil Engineering and Oceanography in oceanographic con- 
struction. The following discussion is intended only to provide some 


background and continuity for the use of the oceanographic data required 











rsa 
to effect a complete solution fo the desten preplem posed. | herrerer. 
ences found in the Bibliography will provide a means for initiating a 


more detailed study concerning spectral analysis of ocean waves. 


Recording Wave Data 


The answer to the rhetorical question of how to extract a 
characteristic frequency from those occurring in a natural phenomenon 
is affected by the methods available for recording the actual physical 
profile of the surface disturbance. The various recording devices 
currently in use all have the common failure of a lack of a standard 
measurement and of distortions inherent in the recording apparatus. 
iMvemertect of the recordinp*distortions becomes more pronounced if 
small non-linear components must be recorded accurately; but the 
presence of these inherent distortions will influence the shape of any 
wave profile recorded and will, therefore, be incorporated in the data 
for a nonlinear analysis as well. Because the recording process is 
quite complicated, the author can only refer to oceanographic literature 
and follow Kinsman (1965) by assuming that for the case in question 
an adequate and representative wave record is available to the 


@esigner from the oceanographer. 
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The Wave Spectrum 


The assumed wave record obtained from the oceanographer 
may also provide an answer to the again rhetorical question of how to 
extrapolate a characteristic design frequency from the natural wave 
phenomenon. The continuous distribution of the frequencies in a wave 
train suggest that a statistical technique may be employed to extract 
the characteristic frequency. The Sverdrup and Munk Theory, for 
example, describes the sea surface by a single sinusoidal component 
called the "significant wave'' (Kinsman, 1965). The next level of 
complexity involves a Fourier analysis. 

In a Fourier or harmonic analysis, the wave record is mechani- 
cally transformed into a spectrogram (McLellan, 1965). If the charac- 
teristic solutions to Eq. (2.1) are orthogonal with respect to the 
weighting function (which arises from the Sturm- Liouville problem 
generated by the separation of variables) and are piecewise differenti- 
able in the recording interval T (Hildebrand, 1962), the solutions 


may be represented by a product series of the form 
7) (t) -) n. (t)- (3.1) 
n 


The solution determined in Chapter 2 was a cosine function. Since 
this function is symmetrical with respect to the z-axis, it may be 


expressed ina Fourier series for an even function; 1i.e., 
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H 
nm {t) = ns = cos ( 0. tiaras eo Jug (Sea) 


n= 1 





Thecoefficients» His anesdetermimedsebyamai le ieee c hed sotaiag. 
n 


(3.2) by the weighting function and by the th 


solution, Be given 
invEq. 3.1, and thenmintegmating over the truncatedsrecoraimesinte ns. a. 


Ae 


n 


00 
H = = n(t) ° cos eee Seal (S05) 


The value given by Eq. (3.3) is equal to twice the average 

value of the product n(t)° cos ( a: - € n) over the recording interval 
T. The wave energy spectrum may also be derived from Eq. (3. 2) 
by squaring both sides of the equation and integrating over the recording 
interval T. The harmonic analysis of the wave record does not pro- 
vide any information regarding the shape of the surface disturbance 
when 

| t | > T (Kinsman, 1965). (3. 4) 
The next logical progression from this discrete spectrum is a repre- 


sentation of the wave record by a continuous spectrum. 
The Bretschneider Continuous Spectrum 


The purpose of this discussion on the methods of presentation 
of wave data is to enable the designer either to request wave data in 


a particular form or to transform a wave record available from the 
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oceanographer. One such method of wave data presentation currently 
in use is the Bretschneider Power Spectrum. This spectrum requires 
that the average height of the highest one-third waves, Elwin as 
determined from a Rayleigh probability distribution of wave heights, 
be known (Dean, 1966). From the value of the significant wave height, 
B73" a parameter which is equal to twice the total energy contained in 
the waves is defined by 


Hy /3 a Deuce NE, . (3.5) 


Bretschneider determined that the spectral distribution of the wave 


energy, S472) as a function of frequency 0 could be represented by 
0 
Syed 0 o 4 
r it Se Oe 5) 
pac (0) = T ( - )" e 0 (2. 6) 


The parameter 0. is the frequency at which the peak energy occurs 


and is defined in terms of the significant wave period, Layee by 


Cc =. (B27) 


Integrating Eq. (3.6) from zero to infinity yields twice the energy 


contained by the waves Ee 
A Functional Design 


The design of the resonant breakwater considered in this thesis 


is an example of a functional design discussed by Wiegel (1964); 








ou 
however, only part of the protection function (which is but one aspect 
of the over-all function of the breakwater) is developed. The functional 
criteria also affect the method required for presenting the wave data. 
The purpose for which the harbor is being comer ructeaaretaics Cour tac 
functional requirements of the resonator and the required method of 
wave data presentation. 

As discussed previously, the wave height and wave period 
parameters normally available from the oceanographer are those of 
the significant wave, an0/3 and a3, The d@stemierltemias wt 
normally require protection from the effects of the highest waves 
Peecent at the structure. The significant wave height, ee has 
been determined to be an accurate index to struetural damage by 
waves (Wiegel, 1964). 

Harbor oscillations or seiching will also normally be found 
in the design criteria for a harbor. Seiching is a complicated response 
of the harbor to a periodic force. This response is a function both 
Smtne geometry of the harbor and of the period of the force. ~ Wave 
data for designing protection from seiching, therefore, must be a 
function of frequency. 

One solution to the problem posed at the beginning of this 
chapter of how to determine a single design frequency Pei the con- 
tinuous spectrum of frequencies occurring in the natural phenomenon 


appears to lie with the statistical significant wave period, We 
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The Bretschneider Power Spectrum derived from the significant wave 
hevohnt fay will then provide a means of evaluating the effectiveness 
of the resonator to an input composed of an infinite and continuous 
spectrum of frequencies. By means of this power spectrum, the effect 
of a resonator on a wave system occurring in nature may now be 


evaluated. 








4. AN EXAMPLE OF A BUNCTIONSA TE Pe rg) 


The Bretschneider Power Spectrum offers a solution to the 
design problems posed at the beginning of Chapter 3. Evaluating the 
effectiveness of the resonator over a spectrum of frequencies by com- 
bining this power spectrum with the reflection and transmission coef- 
ficients represents a linear analysis. The analysis is linear since 
both the power spectrum and the expressions for the coefficients were 


derived from linear solutions. 


A Linear Analysis 


In a linear process, the output, 0, is a linear function of the 
paerut, IL; i.e. , 
Oe— sR (Cee ee (4. 1) 
R(0) is the linear transfer function and is a function of frequency, O. 
The essential features of a linear process are the following (Dean, 1966): 
le Changing the amplitude of the input function, 
I, changes the amplitude ef the eutput function, 
OG proper onally. 
2. The constant of proportionality is the linear 
transfer function and is a function of frequency, 
0, only (i.e., the transfer function is inde- 
pendent of the amplitude). 


38 Input functions may be superimposed. 


The effectiveness of the resonator at a given frequency, therefore, 
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may be determined from the product of the squavecscot 2 a-5( 2.0) son 


(2.19) and the value of the ordinate S 5, from Eq. (3.6) for the 
H 


Mi 
same corresponding frequency; 1.e., 
Se ee (OS 2 (4. 2) 
aa i' oR ay een 
Mr Mi 
or 
Sb .(o.)- a OS ae (4. 3) 
H2 i! a ae a ih 
Mt Mi 


Wave Data anda Linear Analysis 


As an illustrative example of a linear analysis, it is assumed 
that a previously undeveloped coastal area has been selected for the 
construction of a new harbor facility. The harbor has been designed 
to accomodate largeseommercial vessels (lennes and Ekse, 1955) @ The 
proposed entrance to the harbor between the breakwaters is to have 
a depth of fifty (50) feet and a width at the mean-water level of two- 
hundred (200) feet (Figure 4.1). The effectiveness of the construction 
of a single resonator on the local wave conditions is to be measured. 

The preliminary engineering report includes the collection of 
extensive oceanographic data on the proposed location. From this 
report, the significant wave height and period were determined to be 


the following: 
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ae = 12 feet 


"EB = | 257 sec 


The total energy contained in the waves is determined from Eq. (3. 5) 
LoOnDe 
E. = 18 sq. ft. 

The significant wave period, 2) /3? determines the value of os in 
Eq. (3.7) from which the spectral distribution of the incident wave 
enerey, go . may be determined by means of Eq. (3. 6) (Figures 

Ht Mi 
4.2 and 4, 3). 

By satisfying the requirement of Eq. (1.6), the shallow water 
approximation for the wave celerity is determined to be forty (40) feet 
per second through the breakwater entrance. By means of Eq. (2. 2) 
and the value of the wave celerity, the siguneeen wave length, Ls is 
found to be six-hundred (600) feet. Utilizing the design parameters 
from Eq. s (2. 20) and (2. 22), the dimensions of the resonator are as 
follows: 

Lp = 150 feet 

DR = 300 feet 
The values of the reflection and transmission coefficients may now 
be determined as a function of frequency from Eqs. (2.18) and (2. 19), 


respectively. The frequency dependence of the square of these 


coefficients is also shown in Figures 4. 2 and 4. 3. 
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Figure 42. Determination of the effective frequency band-width of the resonator. 
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Figure 4.3. Determination of the effective frequency band-width of the resonator. 
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The reflective and transmissive effectiveness of the resonator 
at a given frequency is determined from the product of the ordinate of 
the incident wave spectrum and the ordinate of the transfer coefficient. 
The resulting spectral distributions of these two linear analyses are 


SO. MOmenloumese te 2 amc .. 
Discussion 


The resonator appears to have a frequency band-width ranging 
from 13 to 17 seconds in which reflection is approximately 100%. The 
effectiveness of the resonator is determined from the amount of energy 
which it reflects and transmits. This amount of energy is measured 
from the area under the power spectrum curves for the reflected and 


transmitted waves, S 93 and § 2 , shown in Figures 4. 2 and 
H Mr H mt 


ao respectively. The narrow spike occurring between the periods 
of 13 and 14 seconds prevents the resonator from being greater than 
85% effective in reflecting the incident wave energy. The proposed 
design for the given wave conditions is approximately 75% effective 


for reflecting the energy contained in the incident: wave spectrum 





between the periods from: 10 to 17 seconds. 
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